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ABSTRACT 



A proof-of-concept experimental validation of a proposed idea for a prototype optical 
system was conducted. This system will be incorporated into a new type of infrared, 
optically multiplexed imaging and multispectral imaging system. This system will use two- 
sided, transmitting-reflecting encoding Walsh masks to form a two-dimensional optical 
Kronecker product. First, a ray tracing design was made to model the optical system. Then 
the optical system was prototyped and ronchigrams were photographed to document the 
aberrations present in the optical system. It was shown that spherical mirrors could be used 
to accurately reimage an object onto the encoding masks without significantly affecting the 
optical accuracy of the image. The geometric aberrations resulting from this design did not 
significantly effect the overall ability to produce the Kronecker product. 
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I. INTRODUCTION 



A. MOTIVATION 

This thesis project is one part of a larger research and development 
program to design and to build a new type of infrared, optically multiplexed 
imaging and multispectral imaging device. This device was conceived by the 
author’s advisor, D. S. Davis. This unique state-of-the-art imaging system will 
use spatially orthogonal Walsh functions to encode the entire image at once, 
rather than use the traditional pixel-by-pixel raster scan method. This new 
approach will allow both spatial imaging and spectral mapping to be 
accomplished with a single device. Thus, Davis’ system will be highly efficient 
and versatile enough to provide imaging and target discrimination tasks in 
infrared remote sensing, astronomy, and surveillance [Ref. 1]. The validity of 
this technique was demonstrated by R.H. McKenzie, Capt., USMC, in an earlier 
thesis [Ref. 2], and subsequent thesis research was conducted by B. J. 
Musselman, LT, USCG, [Ref. 3]; J. P. Sargent, LT, USN, [Ref. 4]; and G. R. 
Parriott, LT, USN, [Ref. 5]. 

B. THESIS GOAL 

Davis’ technique for multiplexed imaging attempts to improve the 
efficiency of previous imaging systems. As described in Reference 1, the new 
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device will encode images using spatially orthogonal masks and simple discrete 
infrared detectors. 

The approach to image multiplexing using orthogonal masks is certainly 
not new [Ref 6, 7]. However, previous image multiplexing devices have two 
characteristics in common: they have employed transmitting-opaque spatial 
masks, and they have required that these masks be manipulated precisely in 
two spatial dimensions. These features represent inefficient, cumbersome 
liabilities in the traditional approaches to multiplexed imaging. The 
transmitting-opaque mask idea necessarily wastes light that strikes the opaque 
regions of the masks, with typical losses approximately equal to 50%. The two- 
dimensional mask manipulation problem has dictated that, to encode n image 
pixels, 2n 2 masks must be precisely positioned by an elaborate 2-D opto- 
mechanical servo system. 

Davis' new approach, as depicted in Figure 1, circumvents both difficulties. 
First, it employs transmitting-reflecting masks and two optical detectors. One 
detector senses reflected radiation, while the other intercepts the transmitted 
portion. Second, it was realized by Davis that it is not necessary to 
manufacture and manipulate a large basis set of masks to multiplex an image 
[Ref. 1]. A mathematically complete two-dimensional encoding basis set can be 
generated by manipulating small basis sets of one-dimensional masks, in 
transmitting-reflecting pairs, and then multiplexing various combinations of the 
encoded signals. 
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Figure 1 Schematic of New Instrument Configuration 



The actual mathematical operation that synthesizes a complete two- 
dimensional mask encoding basis from one-dimensional components is a 
Kronecker product [Ref. 6, 8]. The key to understanding the new multiplexing 
technique is to realize that such a Kronecker product may be performed by 
purely optical means, as shown in Figure 1 and described in Reference 1. All 
that is required is a precise reimaging of one transmitting-reflecting one- 
dimensional mask onto another, from both sides, as shown. When an image 
field of interest is focused on the first mask, the optical scheme does the rest, 
as each 2-D encoding function is generated by cycling through various 
combinations of 1-D masks. There are two potential technical engineering 
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problems that could prevent this scheme from working. The first obstacle is the 
development of a precise one-dimensional mask positioning servo system. This 
problem was solved in the aforementioned thesis research of LT Sargent [Ref. 
4]. The second potential problem deals with the inherent optical accuracy with 
which the two-sided masks reimaging can be achieved. This second problem will 
be addressed by this thesis. 

In an earlier thesis, LT Musselman [Ref. 3] demonstrated that physical 
optics effects (i.e. diffraction) should not have a significant effect on beam 
quality when images are encoded using the new technique. Therefore, it follows 
that successful realization of a working optical Kronecker product instrument 
will depend upon the geometrical optics accuracy with which the two-sided 
reimaging can be done. Therefore, this thesis research has been concerned with 
the following questions: (1) Can precise reimaging (and therefore an accurate 
optical Kronecker product) be accomplished using relatively inexpensive, all- 
reflecting spherical and planar mirrors? (2) How severe are the geometric 
aberrations that result from such a design? (3) How will these aberrations 
affect the overall ability to generate optical products? 



4 



II. Ray Tracing Theory 



The thesis research of B.J. Musselman [Ref. 3] indicated that physical 
optics effects, such as diffraction, will have a negligible effect on the design of 
a prototype multiplexing imaging instrument. Therefore, this thesis project has 
been conducted entirely within the domain of geometric optics. Within this 
realm, light is assumed to follow simple mechanical trajectories, called rays, 
with its propagation governed solely by the basic laws of reflection and 
refraction. 

The advent of the digital computer has lead to the development of efficient 
numerical procedures for calculating the trajectories of rays through complex 
optical systems. The interested reader is referred to other sources for specific 
details on this general theory [Ref. 9, 10]. The commercially available ray 
tracing software used in this thesis is described more fully in Chapter III. 
However, because the optical design of the multiplexed imaging system makes 
use only of simple plane and concave spherical mirrors, it is not necessary to 
delve into the complete theory of ray tracing. Rather, we have developed a 
much simpler approach, which is described below. 

The reader is directed to Figure 2. For the case of a spherical, concave 
mirror, an arbitrary ray originates at the object position (X^ Y 0 , Z G ). We wish 
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Figure 2 Reflection From a Spherical, Concave Mirror 
to calculate the path that the ray will follow when it reflects off the mirror at 
point (X,n, Y m , Z m ) and then intersects a surface (X*, Y it Z s ), as shown in Figure 
2. To simplify matters, the axis of symmetry of the mirror lies on the z-axis and 
the center of curvature for the mirror is at the origin. R is the radius of 
curvature of the spherical surface. All points on the mirror must satisfy the 
condition R 2 = X^ 2 + Y m 2 + Z m 2 . 

The law of reflection governs the ray reflected from the mirror. It is 
expressed in vector notation as, 




\V.\ 1*1 IP, I 1*1 



or, alternatively, 



( 1 ) 



6 



( 2 ) 




where "V* 0 is the vector from (X,,, Y 0 , Z 0 ) to (X^ Y m , Z m ), "V* is the corresponding 
vector from (X*,, Y m , Z m ) to (Xj, Y i} Zj), and Tf is the normal or radial vector from 
the origin to (X*,, Y m , Z J. 

The physical consequences of this law imply two important ideas: 

1. The incident single is equal to the reflected angle ( 0 j = 0 r ). 

2. The normal, incident, and reflected rays are all coplanar. 

Thus, for the cross product in Equation 2 to equal zero, either 0 = 0° or one of 
the two terms equals zero. The first case is trivial and not very interesting. 
However, the second case says that, 



The incident, reflected and normal vectors are expressed in component form as 
follows: 




(3) 



K -(x M -xjf + (r m -rjf + (z m -zj£ 



(4) 




(5) 



N = X m i + Y m j + Z m K 
If the variables A 0 and A; are defined by, 



( 6 ) 
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A . -- ^x m -xf + (T.-rf + K-zf 
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A , = J(X m -X) 2 + {Y n -Yf + (Z m -Zf 


(8) 



this leads to, 
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From Equation 3, the TT is parallel to the left hand side of Equation 9. 
Therefore, each component of TTis proportional to the corresponding component 
of Equation 9, 



X-X, X m’ X i 

—2 — £ + — — i = BX 



(ii) 



Y--Y a Y-Y, 

mo mt 



A o A i 



BY_ 



( 12 ) 



Z.-Z. + V5 

A . A . 



BZ m , 



(13) 



where B is a positive, and as-yet-to-be-determined, constant. 

Numerical ray tracing is based on Equations 11, 12 and 13. The usual 
approach fixes the position of the object and then allows various rays emanating 
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from the object to reflect off the mirror and to terminate on a designated plane. 
The designer formulates the optical system on paper, then inputs the 
parameters into a ray tracing computer program. The program allows several 
of the parameters to be varied in order to reduce the image error, or aberration. 

There are two types of aberrations: chromatic aberrations, which occur 
from the fact that the refractive index (n) is a function of frequency and 
monochromatic aberrations, which arise from imperfections of the mirror and 
from objects off the optical axis (the line connecting the radius of curvature and 
the center of the mirror). The latter of these is the one of concern for this 
project, because there are no refractive surfaces in the design. The 
monochromatic aberrations, to a first approximation, may be further subdivided 
into spherical aberration, coma, astigmatism, field curvature, and distortion 
[Ref. 11]. 

To illustrate the basic principles of ray tracing, two cases will be 
considered. 

1. Assume that a point object is located at the mirror’s center of curvature, 
which is at the origin, X,, = Y 0 = Z Q = 0. Combining this with Equations 11, 12, 
and 13 yields, 




= BX n , ( 14 ) 



and similarly for the Y and Z components. This can be reduced to, 
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( 15 ) 



R 






= BX m , etc. . 



Now let a ray pass through the origin and strike anywhere on the surface of the 
mirror. Obviously, this ray will travel along the normal to the mirror. By the 
law of reflection, it will then pass through the origin again. We are free, 
therefore, to define an image plane that contains the origin, Xi = Yj = Z x = 0. 
Therefore Equations 15 become, 



which specifies the constant B for the system. 

2. As a second example, assume the point object is somewhere on the Z- 
axis, X,, = Y 0 = 0. Equations 11, 12, and 13 then become, 



and similarly for the Y and Z components. From symmetry, the image is 
expected to lie on the Z axis, = Y } = 0). Equations 18 then become, 




( 16 ) 



or 




( 17 ) 




+ 




( 18 ) 
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( 19 ) 



JxK+(Z m -Zf 

The point at which the mirror surface intersects the Z-axis is called the vertex. 
The vertex-object distance is S 0 = Z m - Z 0 , and the vertex-image distance is 
Sj = Z m - Zj. Then Equations 19 reduce to, 



pl« fit* 



2 

R 



( 20 ) 



In the paraxial approximation, where X m 2 + Y m 2 < < S 0 2 and X^, 2 + Y m 2 < <Sj 2 , 



1 1 2 




( 21 ) 



or, 






1 




2 

R ’ 



( 22 ) 



which is the Gaussian mirror equation [Ref 11]. 

When the paraxial approximation does not apply, different circular areas 
at the mirror, corresponding to + Y m 2 = constant = r m 2 , cause objects to be 
focused at differing values of S;: 
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Expanding this equation using the binomial series, 



( 23 ) 
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For an object at oo, 





\ 






2 

R ' 





2 

R ’ 



or, 
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(24) 



(25) 



(26) 



which is the third-order spherical aberration [Ref. 11]. The additional lo west- 
order aberrations may be derived by similar methods. 

Historically, it was necessary to perform such calculations for each 
aberration by hand, and then to recombine them ray-by-ray in order to discern 
their effects on the overall system performance. Modern ray tracing has 
removed this chore. We now simply program the computer to predict 
performance, with all aberrations folded into the calculations at each step. This 
is the procedure followed throughout this thesis. 
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III. OPTICAL DESIGN PROCEDURES 



A. INTRODUCTION 

In the previous chapter, the basic geometric calculations needed to predict 
ray propagation were discussed. The present chapter discusses how these basic 
procedures are combined to generate an optimized optical imaging system. The 
optimization carried out as part of this thesis research are by no means the best 
possible characterizations of a reimaging system; rather they represent the best 
design for a system using readily available, off-the-shelf spherical mirrors. It 
is virtually certain that better performance could be expected from a system 
employing customized aspheric optics, but such a system would be much more 
costly then present financial constraints would allow. 

The design optimization procedure is straight forward. A geometric object 
is specified, in general terms, and numerical ray tracing is performed by a 
computer. Each object point emits rays that pass through various parts of the 
system. An image plane is specified at the system’s output, and the program 
yields the Cartesian coordinates at which each exiting ray intersects this image 
plane. Ideally, a perfect imaging system would reimage all rays originating from 
a given object point into an infinitesimal image plane. However, real optical 
systems, and the computer models that represent them, exhibit aberrations 
which cause the exiting rays to intersect the image plane in the neighborhood 
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of the ideal point image, rather than at the point itself. When the intersections 
of many such rays are plotted, the density of intersection points specifies a point 
spread function in the image plane. The smaller the point spread function , the 
better the system approximates the ideal. Referring to Figure 1, the reader will 
note that the multiplexed imaging system design requires two basic reimaging 
operations. First, an initial input object (the "chopper" in Figure 1) must be 
imaged onto the first encoding mask. That mask is then reimaged onto a 
second encoding mask, to form a Kronecker product [Ref. 6,8]. The 
optimization of this design requires that the basic characteristics be specified, 
such as radius of curvature and aperture, for the spherical mirrors used. 
Hypothetical object points are then placed in the input plane, and numerical ray 
tracing generates the point spread functions in the designated output plane. 
The program varies the location(s) and orientation(s) of the mirror(s) in order 
to minimize the overall point spread function. It then gives the user some 
numerical parameter, such as RMS ray scatter about a mean image location, 
that characterizes the point spread function and, hence, predicts image quality. 

The commercially available ray tracing software used in this thesis has 
been the BEAM FOUR system, distributed by Stellar Software of Berkeley, 
California. The following sections detail the use of BEAM FOUR in the 
optimization of the multiplexed imaging instrument’s optical design. 
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B. BEAM FOUR 



The first phase of this thesis research consisted of designing the optical 
system to image an object through the two sets of transmit ting-reflecting masks, 
as outlined in the previous section. BEAM FOUR is an IBM PC compatible 
program that performs ray tracing in Cartesian space. It requires both an 
optics table and a ray table to generate a ray trace with twelve significant 
figures of precision. Optics tables specify the optical elements to be used and 
ray tables indicate how the optical system is to be illuminated. BEAM FOUR 
has an AUTOADJUST command that will optimize the focus of the system by 
altering designated parameters. Furthermore, the program is interactive. The 
optical system, and rays, can be viewed with the LAYOUT command. The 
(X,Y,Z) coordinates at which rays intersect optical surfaces can be listed with 
the VERTEX command. There is also a the PLOT command, which is useful 
in creating spot diagrams to help analyze the optical system’s performance by 
giving a visual representation of the point spread function. 

Optics tables are a list of the optical elements to be analyzed by BEAM 
FOUR. Each optical element is listed in the sequence that it will be 
encountered by the light. All the characteristics (position, orientation, radius 
of curvature, etc.) of each optical element is specified in the table under the 
appropriate header. All optics disc files table are identified with file name 
extensions of ".OPT". The file preamble lists the number of optical surfaces and 
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the table title on the first line; each column’s content on the second line. The 



following is a list of column headers used: 



Index: 

Diameter: 

X,Y,Z: 

Tilt: 

Pitch: 

Roll: 

Curvature: 



Figure: 

Mirror: 



index of refraction of the medium in which rays approach 
a surface 

outer diameter or aperture of the optical component 
the vertex coordinates of surface 

degrees of rotation of the optical element around the X axis 
degrees of rotation of the optical element around the Y axis 
degrees of rotation of the optical element around the Z axis 
Value = l/(radius of curvature); positive means curved 
towards the +Z axis, negative means curved towards the 
-Z axis, zero means a flat surface 
round or square geometrical form 
surface identity (mirror, lens, iris) 



If one of the optical element’s characteristics is not designated, it will use a 



built-in default value. The default value for the index of refraction is 1.00 and 



for all other numbers is 0.00. The Z axis is the optical axis. 

Ray tables are set up similarly. Each ray’s position of origin, initial 
direction, goal position and final position are designated under the appropriate 
header. Ray table disc files are identified with file name extensions of ".RAY'. 
The preamble lists the number of rays and the title on the first line; each 
column’s content on the second line. The following is a list of column headers 
used: 



X0,Y0,Z0: 

U0,V0,W0: 

Xg,Yg,Zg: 

Xf,Yf,Zf: 

@wave: 

Note: 



coordinates of ray starting point 

direction cosines at which ray vectors leave X,,, Y 0 , Z 0 

X,Y,Z goal coordinates at the final surface 

X,Y,Z actual coordinates at which the ray intersects the 

final surface 

changes color of the ray 

indicates the surface at which the ray terminated 
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Coordinate default values are 0.0, and initial ray directions that are not 
specified propagate parallel to the Z-axis in the positive Z direction. 

The AUTOADJUST optimization feature uses a Gauss-Newton nonlinear 
iterative least-squares routine to automatically adjust tagged parameters for the 
least-mean-square discrepancy between the final ray state and user-specified 
goal. This process works iteratively, and exits when successive reductions in 
error are smaller than e* 8 . Up to ten independent parameters may be optimized, 
the most common being surface location, pitch and radius of curvature. The 
root- mean- square (RMS) error in the fit is displayed for each iteration so that 
the operation’s progress can be monitored. 

A quick and easy way to verify that the trace has done the intended job is 
to examine it pictorially. The LAYOUT command creates a three-dimensional 
on-screen diagram of the optical system, complete with rays. The diagram can 
be rotated about any axis, and there is a zoom feature that allows details to be 
viewed. 

A ray vertex is the point where a ray segment joins the following ray 
segment by being reflected or refracted at an optical interface. The LAYOUT 
command gives a quick view of the approximate ray positions, but for a more 
accurate position the VERTEX LIST command is used. It lists the (X,Y,Z) 
coordinates of each ray at each optical surface. 

The PLOT command displays the spatial relationship of the rays at the 
final surface on a two-dimensional plot. This helps reveal first-order errors in 
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focusing. Once the system has been refined, a spot diagram can be created 
using the sub-command RANDOM which uses a Monte Carlo random ray 
generator to fill the plot. 

Ce CONSTRAINTS 

In designing the optical system, the following constraints were used in the 
design process: 

1. The final multiplexed imaging device will produce images using radiation 
with radiation up to 20 pm in wavelength. Taking diffraction effects into 
account, the RMS point spread deviation for the system should not be more 
than approximately 2 - 3 times the longest wavelength (i.e. about 50 pm, 
corresponding to BEAM FOUR parameter RMS dev = 0.0050). 

2. Since the device will be used in remote sensing applications to view 
objects at long distances, assume that the maximum angle off the optical axis 
at which rays will be propagated is ±0.5 arcminutes. 

3. When completed, the imaging device will have applications for 
shipboard, aircraft, or satellite use. Therefore, its physical dimensions should be 
restricted to about 1 cubic meter. 

4. Due to budget considerations, the cost of the optical elements should be 
minimized. Spherical mirrors are readily available through commercial sources. 
Therefore, their use is preferred over other shapes in order to be cost effective. 
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5. All calculations in BEAM FOUR will be done assuming centimeters as 
units. 

D. INITIAL DESIGN 

Figure 1 shows the symmetry about the plane containing the chopper and 
the two masks. Taking advantage of the system’s symmetry, only one side of 
the system needs to be designed then the other side will be the mirror image. 
The initial design configuration is shown in Figure 3. The line drawn between 
the center of two optical elements, usually following the optical axis, is the 
principle ray. Figure 3 has six optical elements. The principle rays connecting 
these elements will be designated the principle axis. Nine rays of light (one on 
the optical axis, four spread vertically and four spread horizontally) parallel to 
the optical axis were used to simulate a point source at infinity. A simple 
biconvex lens with index of refraction, n = 1.5 is use to focus the light onto the 
injector mirror. The light is then focused onto the first mask with a spherical 
mirror. The optimized solution for this configuration is listed in Tables 1 and 
2 (HUGA.OPT and HUGA.RAY), and has an RMS point spread deviation = 
0.0010 cm. All subsequent BEAM FOUR design configuration tables are 
contained in Appendices A, and will be referred to by their preamble title. 
When the second spherical mirror is added to reimage the point source 
optimally onto the second mask, the system has an RMS point spread deviation 
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Mirror i 



Mirror 2 




Figure 3 Initial Design Configuration 

of 0.0033 cm. This system’s parameters are listed in HUGB.OPT and 
HUGB.RAY. 

The next design step was to look at an object located off the optical axis. 
For the system concerned, locating an object ten meters from the focusing 
element will approximate an image at infinity. Using the assumption that the 
object is a maximum of 0.5 arcminutes off the optical axis translates to a lateral 
displacement of 0.1454 cm at 10 m. Thus, using the previous optical layout in 
HUGB.OPT and offsetting the point source 0.1454 cm above the optical axis 
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TABLE 1 



OPTICS TABLE FOR THE INITIAL DESIGN 



5 surfaces 




HUGA. 


.OPT 








Inde:< Diameter 


X 


Z 


Pitch 


Curvature 


F iqure 


Mx rror 


: 3.0 : 




9.85 : 




0.0074074 : 


s Circle : 


lens 


1.5 : 3.0 : 




9.95 : 




-0.0074074 : 


: Circle : 


lens 


: 0.5 : 




145.0 : 


-45 . 0 




: Square : 


mirror 


: 5.0 s 


36.6091 


145.0 : 


94.6314 


-0.0270400 ! 


: Circle : 


mirror 


: 5.0 : 


-0.17678 


151.00 : 


90 . 0 




: Square : 


mirror 



TABLE 2 

RAY TABLE FOR THE INITIAL DESIGN 

9 rays HUGA.RAY 



XO 


YO 


X f i na 1 


Yf i na 1 


Zf i na 1 


Xqoa 1 


Yqoal Zqoa 1 


Note 


1.2 




-0 . 1 768 


: 0 . 0000 


: 151.0015 


-0. 17678 


: 151.0 


: o! 


5 


0.6 




-0 . 1 768 


: 0 . 0000 


: 151. 0004 


-0 . 1 7678 


: 151.0 


: ok 


5 


0 . 0 




-0. 1768 


: 0 . OOOO 


: 150.9994 


-0.17678 


: 151.0 


: ok 


5 


-0 . 6 : 


! 


-0. 1768 


: 0 . 0000 


: 150.9984 : 


! -0.17678 


: 151.0 


: ok 


5 


-1.2 




. -0. 1768 


: 0 . 0000 


! 150.9972 : 


i -0.17678 


: 151.0 


: ok 


5 




1 . 2 


-0 . 1 768 


: -0.0022 


: 150.9994 : 


! -0.17678 


: 151.0 


: ok 


5 




0.6 


-0. 1768 


: -0.0011 


: 150.9994 ; 


: -0.17678 


: 151.0 


: ok 


5 




-0.6 ! 


! -0.1768 


: 0.0011 


: 150.9994 : 


: -0.17678 


: 151.0 


:ok 


5 




-1.2 ! 


: -0 . 1 768 


: 0 . 0022 


: 150.9994 : 


: -0.17678 


: 151.0 


:ok 


5 



yielded an RMS point spread deviation = 0.5813 cm (HUGC.OPT and 
HUGO. RAY). This was unsatisfactory, being over hundred times the allowable 
error specified above. 



E. ALTERNATE TECHNIQUES 

Next, one or both of the spherical mirrors was replaced with aspherical 
mirrors having more exotic shapes using the SHAPE command in BEAM 
FOUR. Their positions were also allowed to vary. This yielded only minor 
improvements in reducing the overall point spread error. For object points on 
the optical axis, the error was reduced by only 1-3 pm. When the source was 
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off the optical axis, the RMS point spread deviation was typically 0.0560 cm. 
This corresponds to a reduction in error by a factor of ten from the initial trials, 
but it is still too large. 

Davis speculated that the major source of error was in the external 
biconvex lens focusing element used to simulate a telescope. He suggested using 
a Petzval lens system to reduce the error when focusing the light. Using the 
Petzval lens system [Ref. 12:p. 429] to replace the biconvex lens yielded 
improvements of the same order of magnitude as those achieved when using 
aspherical mirror surfaces. This degree of image point spread was still too 
large. 

F. FINAL SOLUTION 

The persistently large RMS beam spread error present when imaging a 
point off the optical axis further suggested that this error was caused by the 
inclusion of the external telescope focusing element. Therefore, in order to test 
this hypothesis, the external element was eliminated from the design. This is 
practical since Davis’ imaging system itself does not included the external 
telescope. In eliminating the telescope, the assumption that the radiation 
propagates through the system in such a way that the object is imaged perfectly 
at the chopper was made. 

Figure 4 shows a diagram of the plane containing the chopper and the two 
disks with the Walsh encoding masks. Mirror 2 images the first mask onto the 
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second mask, and this mirror’s optical axis must bisect the angle formed by the 
two masks and mirror 2. This requirement arises from the fact that the masks 
are transmitting and reflecting, in order to generate an optical Kronecker 
product. Therefore, all the radiation must be reimaged on the second mask. 
Hence, the first step is to find the optimum position of mirror 2 that gives an 
acceptable beam spread error at the second mask. By fixing the position of the 
object (point source) and the image at the centers of mask 1 and 2 respectively, 
and letting the position of the mirror 2 and its radius of curvature vary, the 
results in Table 3 were calculated. The radii of curvature (R) listed correspond 
to standard mirror values that are readily available commercially. 

Using the mirror with R = 50 cm, if the point source is moved to the edge 
of the mask (assume the mask is 5 mm 2 ), the RMS point spread deviation 
becomes 0.0025 cm (HUGD.OPT and HUGD1.RAY), which is within the desired 
limit. This confirmed that the original error was caused by the external 
telescope focusing elements. Hence, those elements were permanently excluded 
from the design. 
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CHOPPER 




Figure 4 Plane Containing the Chopper and Two Encoding Masks 
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TABLE 3 



RMS POINT SPREAD DEVIATION FOR SEVERAL MIRRORS 



Radius of Curvature 


RMS Deviation 


20 cm 


0.0054 cm 


30 cm 


0.0035 cm 


40 cm 


0.0026 cm 


50 cm 


0.0021 cm 


60 cm 


0.0017 cm 



Next, the position and radius of the mirror 1 must be determined. To 
avoid having mirror 1 and 2 physically too close, which might cause mechanical 
problems when the system is assembled, mirror l’s radius of curvature was 
chosen to be 60 cm. Analysis of the geometry of the system indicates that this 
mirror must lie on the extended principle line between mirror 2 and mask 1. 
This principle line was made the optical axis (Z axis) in the BEAM FOUR by 
rotating the plane containing the chopper and two encoding masks 2.8° from 
vertical. The coordinates of the centers of the chopper input, mask 1 and 2 are 
known from Figure 4. When all this information was input to BEAM FOUR an 
optimized position of mirror 1 was calculated (HUGE. OPT and HUGE. RAY). 
An RMS point spread deviation of 0.0041 cm was achieved at mask 1. When 
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this is combined with the position calculated for mirror 2, an overall system 
RMS point spread deviation = 0.0023 cm (HUGF.OPT and HUGF.RAY) for the 
image at mask 2 was calculated. When the object was moved to the edge of the 
mask, the total system RMS point spread deviation was 0.0042 cm (HUGF.OPT 
and HUGF1.RAY). Thus, this system meets the requirements for the maximum 
allowable point spread error cited previously in section C of this chapter. 
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IV. PROTOTYPE OF THE OPTICAL SYSTEM 



A. RONCHI TEST PROCEDURES 

Experimental evaluation of the actual prototype design required that some 
objective procedures be adopted for judging image quality. Ideally, a precision 
interferometric and/or Fourier optics analysis of system performance would have 
been employed. However, the elaborate apparatus needed to do such analysis 
was not readily available. Therefore, a conventional Ronchi ruling test was 
chosen. 

Until recently, there were historically two general ways of thinking about 
image quality in optical systems. The first of these was a purely geometrical 
optics approach, and was based upon analysis of ray shadow patterns, such as 
Focault test, Ronchi test, and their relatives [Ref. 13]. The second approach 
was purely in the domain of physical optics, and involved interferometric 
analysis of wavefronts exiting from the optical system under test [Ref. 11]. The 
recent advent of Fourier optics [Ref. 14] has lead to an elegant melding of the 
two historical approaches. The approach used for this project employed a simple 
form of the Fourier optics point of view in our system analysis. 

To understand this approach, consider first the properties of an ideal 
imaging system, as it images a point source (a spatial Dirac delta function). If 
diffraction were not present, the ideal image would be another spatial Dirac 



27 



delta function. However, because of finite apertures and optical wavelengths, 
what actually appears in the image plane is a blurred distribution of irradiance, 
E(x,y). This distribution is called the point spread function, because it describes 
the image spreading characteristics produced by an isolated point source. 

Basic diffraction theory [Ref. 11] shows that the actual distribution of 
irradiance within the ideal system’s point spread function is related to the two- 
dimensional spatial Fourier transform of the limiting aperture stop, or pupil, of 
the system. Hence, a system with narrow apertures will exhibit a wider point 
spread function than will a system with a wide apertures. For systems with a 
rotational axis of symmetry, it is a common practice to simplify the analysis to 
a one-dimensional cross-section of the point spread function along a diameter, 
and to threat this as a corresponding one-dimensional radial Fourier-Bessel 
transform of the limiting aperture. This cross-section is called the line spread 
function. 

When an extended object, consisting of a continuum of continuous points, 
the situation becomes slightly more complicated. If the point spread function 
is invariant across the field of view containing the image, then each image point 
will be broadened by the same point spread function. The observed image 
irradiance distribution will be the weighted sum of all these point-by-point 
contributions. It can be shown [Ref. 11] that this irradiance distribution is the 
convolution of the point spread function with the ideal distribution that would 
result in a diffraction-free (geometrical optics limiting) case. The same 
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conclusion holds for the one-dimensional line spread function concept. It is, 
therefore, apparent that a convenient way to quantify the quality of an image 
produced by a system is to measure the point spread and/or line spread 
functions of that system. Since the point spread and line spread functions 
always manifest themselves as components in a spatial convolution, it is also 
common practice to express them in Fourier transform form, so that the 
convolutions may be considered using the Fourier convolution theorem. The 
Fourier transforms are called optical transfer functions, and their modulii are 
the modulation transfer functions. To analyze the prototype system, the direct 
point/line spread functions, rather than the transfer functions were used. A 
reasonable measure of image quality is the RMS distribution, or standard 
deviation, of irradiances about the point or line spread function’s mean 
coordinates. 

When optical aberrations are present in the system, the point and line 
spread function ideas may be extended to include those effects. Aberrations 
will, in general, introduce asymmetries and distortions into the spread 
functions. They may also destroy the translational invariance of the functions, 
but this problem may be minimized by considering only narrow fields of view, 
as was the case in this thesis research. In this event, the line and point spread 
functions are still spatial Fourier transforms, but they are transforms of virtual 
apertures, with complex spatial characteristics. In any event, the resulting 
spread functions are still useful measures of image quality. 
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For this thesis research, we have chosen to quantify system performance 
using the RMS width of the line spread function. This proved to be convenient, 
because the line spread function can be obtained relatively easily from a 
conventional Ronchi test. 

Figure 5 shows a typical arrangement used to perform a Ronchi test. The 
centers of the light, the mirror and the reticle must all be in the same plane. 
The slit is placed in the plane containing the reticle, with the slit parallel to the 
reticle. When the slit is illuminated, a concave mirror will reflect the light and 
form an image of the slit on the reticle. If the slit width is substantially 
narrower than the spacing of the reticle rulings and the mirror is placed at a 
distance equal to its radius of curvature from both the reticle and the slit, the 
resulting image magnification is minus one. Thus, when the image is centered 
on either a dark or clear reticle line, the system’s exit pupil (the mirror) seen 
when looking through the reticle will be uniformly dark or light, if the line 
spread function is substantially narrower than one half a reticle ruling spacing. 
When the mirror is at a position other than the radius of curvature, the 
magnitude of magnification for the image formed will not be one. The resulting 
image’s line spread can be greater than the width of the reticle rulings, and the 
exit pupil will appear to be crossed by light and dark contours, or fringes. Thus, 
the Ronchi test is a sensitive procedure for focusing the system (by minimizing 
the line spread function) and for establishing an experimental upper limit of the 
line spread function. [Ref 11] 



30 




Figure 5 Setup for the Ronchi Test 

B. PROCUREMENT 

Based on the final solution for the basic optical design problem as 
described in Chapter III. F., the prototype for the optical system will use two 
spherical mirrors with radii of curvature R x = 60 cm and It* = 50 cm and 
diameters D x = 5.04 cm and D 2 = 2.54 cm. However, when an attempt to 
procure the mirrors was made, it was discovered that none of the major 
companies had these mirrors in stock. Therefore, due to the limited time the 
author had to complete this project, two spherical mirrors with R! = 152 cm (60 
inch) and = 127 cm (50 inch) were obtained from Edmund Scientific Co. By 
using these mirrors, the same design can be used, but all physical dimensions 
have to be multiplied by 2.54 to convert from centimeters to inches. Another 
adjustment had to made to compensate for these mirrors. Originally, the 
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mirrors in the design were to have diameters, D x = 5.08 cm and D 2 = 2.54 cm, 
corresponding to f-ratios of f/5.9 and f/9.8, respectively. The mirrors available 
had D! = 8.89 cm (f/8.6) and D 2 = 6.35 cm (f/10). Mirror holders had been 
purchased ahead of time to fit the mirrors in the original design. Therefore, 
adapters were made in the physics fabrication shop that support the mirrors 
using the original mirror holders. A drawing of one of the adapters is shown 
in Figure 6. The discrepancies between the actual mirror’s f-ratios and those 
of the original design also required modification of the detailed design 
calculations to predict system performance. 

Use of alternate mirrors caused yet another problem. In anticipation that 
the system would incorporate mirrors with radii of curvature approximately 50 
cm, Davis had purchased commercially manufactured Ronchi rulings with 
density as high as could be accommodated by such a system. However, when 
the mechanical dimensions were scaled up to accommodate the alternate 
mirrors, those Ronchi rulings were too dense. Therefore, two sets of Ronchi 
rulings were generated at densities of one and two line cycles per millimeter. 
The process was a two step operation. First, AutoCAD software was used to 
produce precise tracings on a high definition Hewlett-Packard 7550A plotter. 
Second, the resulting plots were transferred to transparencies, using a 
conventional office copy machine. 
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Figure 6 Mirror Adapter 
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C. DETERMINING THE RADII OF CURVATURE 



As stated in the previous section, the prototype optical system had to 
incorporate off-the-shelf mirrors whose characteristics were quite different from 
the those specified by the original design. Before assembling and testing the 
optical system, the actual radii of curvature for each mirror had to be measured 
and verified because the manufacturer, Edmund Scientific, guaranteed the radii 
of curvature tolerances to only ± 2%. The version of the Ronchi test described 
above was used to make these measurements. 

A one line per millimeter Ronchi reticle was placed at a distance 
approximately equal to the radius of curvature from the mirror. A light source 
with a variable intensity was placed behind the slit to illuminate the mirror 
through the slit. With the slit open several millimeters and high light intensity, 
the mirror was tilted and translated to form an image on the center of the 
reticle. While looking through the reticle at the exit pupil (mirror surface), the 
intensity of the light was reduced and the slit was stopped down (width was 
reduced) until an image with several dark and light lines appeared across the 
exit pupil. This is called a ronchigram. Then, the mirror was mounted on a 
micropositioner so that it could be translated both perpendicularly to the 
reticle’s rulings and, independently , parallel to the reticle normal. Using the 
micropositioner, the mirror was moved until all the fringes disappeared and a 
uniform field of illumination was seen. This was the position where the mirror 
was focused, because the line spread function was sufficiently narrow. The 
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focused position of the mirror was found to be uncertain to ±0.5 mm (given 
estimated uncertainties in the relative positions of the slit and Ronchi reticle). 
This corresponds to a radius of curvature uncertainty of ±0.1 cm. 

Using the one line per millimeter reticle, the ronchigram showed no 
curvature of the fringe pattern in the out-of-focus configuration. With the two 
lines per millimeter reticle, the fringes began to show a small amount of 
curvature, indicating various aberrations in the mirror. 

To summarize, the radii of curvature for both mirrors were measured and 
found to be: 

Rj = 128.7 ± 0.1 cm 
Ra = 152.3 ±0.1 cm. 

D. OPTICAL SYSTEM ALIGNMENT 

Given the actual mirror radii of curvature, it was then necessary to update 
the basic optical system design to more accurately model the laboratory setup. 
Entering the actual radii of curvature for mirror 1 and mirror 2 into BEAM 
FOUR, the optimum coordinates of all the system’s optical element’s centers 
were determined. These are detailed in files HUGG.OPT and HUGG.RAY and 
are summarized in Table 4. Using these values, a working prototype optical 
system was assembled on a conventional flat optical table. The center of mask 
1 is at the origin of the system and is located 11.7 cm above the optical table. 
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The X axis is parallel to the shorter axis of the optical table, the Y axis is 
vertical, and the Z axis is parallel to the longer axis of the optical table. 

Table 4 

COORDINATES FOR THE OPTICAL SYSTEM 





X (cm) 


Y (cm) 


Z (cm) 


Light Source 


7.6 


23.4 


- 0.5 


Mirror 1 


0.0 


0.0 


151.0 


Mask 1 


0.0 


0.0 


0.0 


Mirror 2 


15.2 


0.0 


128.0 ; 


Mask 2 


15.2 


0.0 


- 0.9 



The optical procedures are as follows: 

First, position and align the light source, mirror 1 and mask 1 using simple 
mechanical measuring methods (i.e. a meter stick). Replace mask 1 with the 
Ronchi ruling reticle and perform the Ronchi test to achieve optimum focus (see 
Figure D-l). Because the angle separating the light source and the reticle 
somewhat exceeds the regime where the paraxial ray approximation holds, the 
image quality is poorer than that achieved earlier during the radii of curvature 
measurements. Note that the tilt of the fringes off vertical is approximately 
parallel to a line drawn from the centers of the light source and the reticle. 
This is evidence for astigmatism and coma in the image. Approximately three 
fringes are visible at the position of best focus. 
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Next, replace the Ronchi reticle with a flat mirror and position mirror 2 
and the reticle in place of mask 2. Again perform the Ronchi test and adjust 
the position of mirror 2 until the optimum focus is achieved (see Figure E-l). 
Again, approximately three fringes are visible, implying about the same amount 
of off-axis image aberration as was present in the previous step. 
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V. ANALYSIS OF THE RESULTS 



A. RONCHIGRAM RECORDS OF EXPERIMENTAL RESULTS 

To aid in the analysis of the results, the ronchigrams seen at mask 1 and 
mask 2 were photographed. Davis suggested that the Ronchi ruling reticle and 
slit be rotated in the vertical plane at 30° increments to show how any system 
optical aberrations might vary with angle. These images appear in Appendices 
D and E. Pictures were also taken of the images formed when determining the 
radius of curvature for mirrors 1 and 2, and are contained in Appendices B and 
C. Several pictures were also taken when mirrors 1 and 2 were moved away 
from optimum focus to show how this affected the number of ronchigram 
fringes seen. These pictures were taken with black and white TMAX 400 
KODAK film, with an aperture setting of 5.6, and shutter speed of 1/4 second. 
The camera used was a Nikon F3 with a 80 - 200 zoom lens. 

B. THEORETICAL APPROXIMATIONS OF SYSTEM LINE SPREAD 

FUNCTIONS 

As stated before (see section III. A.), ray spot diagrams can be very useful 
in analyzing optical systems. An example of a typical spot diagram, produced 
using a spherical mirror with a point source and the image plane located slightly 
beyond its center of curvature, is shown in Figure 7. Using BEAM FOUR, spot 
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Figure 7 Spot Diagram 

diagrams were calculated. When many numerically generated rays, uniformly 
distributed in space and angle variables, is allowed to "illuminate" a numerical 
model of an optical system, the local density of spots on a spot diagram should 
be a reasonable predictor of expected radiative flux in that vicinity. By using 
such an approach, the theoretical predictions for the width of the line spread 
functions were calculated that should characterize the actual optical system, 
including all first and higher-order geometrical aberrations. When a large 
number of rays were used, the RMS scatter of spots about the beam centroid 
was chosen as a reasonable measure of predicted line spread. The two lines per 
millimeter Ronchi ruling reticle used for the experiments will produce 



39 



observable fringes whenever the line spread width exceeds 0.025 cm. If the 
RMS line spread deviation for the spot diagrams is assumed to represent the 
line width produced at the designated surface, and proper geometric scaling is 
taken into account, then the number of fringes that should be seen in each 
experimental ronchigram can be predicted, at least approximately, on theoretical 
grounds. 

To get the most accurate prediction, several hundred rays should be used 
to fill the modeled mirror with light. BEAM FOUR is limited to 99 rays per 
calculation, which makes the statistical analysis of the spot diagrams somewhat 
dubious. Therefore, three sets of independent rays were used to search for 
trends and to check that the predicted number of fringes approached the actual 
number of fringes seen in the ronchigrams. 

C. SUMMARY AND COMPARISON OF THEORY AND EXPERIMENT 

As stated in the previous section, three set of rays were generated to 
attempt to develop a trend in the predicted number of fringes seen in the 
ronchigrams. The first sets of rays (P45) consisted of five point sources in a 
line, with nine rays of light (one ray striking the center of the mirror, four rays 
fanned vertically, and four rays fanned horizontally) emanating from each point 
source. The second sets of rays (P25i) used the initial set of rays without the 
four rays from each point source closest to the edge of the mirror. The third set 
of rays (P25o) copied the original set of rays without the four fanned rays from 
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each point source closest to the center. Table 5 contains the results from these 
three simulations. Appendix F contains the BEAM FOUR optical tables and ray 
tables used to produce these line spread predictions. Appendix G contains the 
plots of the line spread scatter diagrams obtained for the ray tables in P45. 

TABLE 5 

PREDICTED AND ACTUAL NUMBER OF HALF FRINGES 





P25o 


P25i 


P45 


Actual # 


Rj Mirror 1 


0.3 


0.2 


0.2 


< 1 


R, Mirror 2 


0.3 


0.1 


0.2 


< 1 


Mask 1 Vertical Slit 


8.3 


7.0 


7.8 


3 


Mask 1 Horizontal Slit 


8.2 


6.8 


7.7 


2 


Mask 2 Vertical Slit 


10.9 


10.3 


8.3 


3 


Mask 2 Horizontal Slit 


4.1 


2.1 


3.4 


2 



Although the ray tables used for P25o and P25i both have the same 
number of rays, the predicted numbers of half fringes for these calculations are 
different. This occurs because the rays striking the mirror closest to its center 
most closely adhere to the paraxial approximation, which minimized aberration- 
induced line spreading. If a ray table with several thousand rays was used to 
calculate the predicted number of half fringes, the predictions should converge 
to the observed number of half fringes. 
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D. SOURCES OF ERROR 



Throughout the prototyping phase of this project there were several 
potential sources of error that had to be considered before the final analysis of 
the system could be completed. The two major sources were the ability to 
position the optical elements accurately, and the optical quality of the reticle 
used in the Ronchi test/procedure. 

1. Positions of Optical Elements 

Although BEAM FOUR was able to calculate the optimum position 
of each optical element to twelve significant figures, in practice, the optical 
elements’ position could only be measured to within one millimeter when they 
were positioned on the optical table. This uncertainty affected the 
measurements taken when the radii of curvature were determined for mirror 
1 and mirror 2. The problem was further compounded when the optical 
elements were positioned, based on the radii of curvature that were determined. 
To achieve the best results, a three-dimensional (X,Y,Z) micro-positioner that 
can also adjust the tilt and pitch of the mirrors should be used. Also, a more 
accurate means of determining the optical elements’ locations relative to the 
designated origin should be used. Only coarse measurements, using a meter 
stick with millimeter rulings, were made for this project. 
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2. Ronchi Reticles 



When the alternate mirrors were substituted for the mirrors that 
were originally specified, the reticles that were originally purchased were no 
longer useful because of their high ruling density. Therefore, "home-made" 
reticles were generated to use in the Ronchi test and positioning of the spherical 
mirrors in the optical system. Due to both the manner in which the ink was 
absorbed onto the paper to create the rulings and to the ability of the 
photocopying process to reproduce the lines on a transparency sheet, the home- 
made reticles were far from perfect. The rulings were not uniform in width, 
and their edges were extremely jagged. These imperfections affected the 
ronchigrams photographed in Appendices B - E. Thus, the widths of observed 
Ronchi fringes must be considered to be less accurate than would be the case 
had high quality reticles been available. 

E. INTERPRETING THE ABERRATIONS 

The ronchigrams produced for this project show several different 
characteristics. A rough, bumpy (see Figure C-l) texture occurs when the 
entrance slit for the light source is very narrow and the camera is focused 
extremely well. This mottled texture can be produced by two sources: (1) 
mirror surface imperfections, and (2) blemishes in the Ronchi reticles. Both 
effects appear here. The half moor! effect (see Figure C-2) is produced by 
astigmatism and coma. The Yin-Yang effect (see Figure D-6) arises because the 
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system uses a three-dimensional off-axis optical configuration in which the 
spherical mirrors "twist" the light about the principle rays. It is strongest near 
the center of the mirrors. [Ref. 13] 

1. Rj Mirror 1 

The uniform dark field seen in Figure B-l indicates that: (1) the 
mirror is at a distance equal to its radius of curvature from the reticle, (2) that 
the line spread function is substantially narrower than one Ronchi ruling cycle, 
and (3) that the image is centered on a dark reticle ruling. Figure G-l shows 
the simulation of the line spread for this situation. The nearly uniform spacing 
of the dots indicate that there should be no major aberrations. Figures B-2 
through Figure B-4 show that moving the mirror ±5.0 mm from optimum focus 
will increase the number of observed Ronchi fringes by one. No significant, 
unexpected aberrations are evident in this mirror; it appears to be a perfect 
sphere to within the tolerances that this test can measure. 

2. R -2 Mirror 2 

The essentially uniform field in Figure C-l indicates that the mirror 
is approximately at a distance equal to its radius of curvature from the reticle, 
as was the situation in the previous case. Figure G-2 shows the simulation of 
the line spread for this situation. The nearly uniform spacing of the dots 
implies that there should be no major aberrations affecting the image. 
However, the rounded fringes in Figure C-2 and Figure C-3 indicate that a 
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slight astigmatism exists in this configuration. It may be that this astigmatism 
is more apparent in mirror 2 than in mirror 1 because of the larger off-axis 
angle between the light source and the reticle (1.5° and 1.2° respectively) for 
mirror 2. 

3. Partial System Containing Mirror 1 and Mask 1 

Figures D-l through D-6 show how the aberrations of the optical 
system affect the light propagating through the system with the entrance slit 
and reticle at various orientations about the principle ray. 
a. Figure D-l 

The entrance slit and the reticle were both vertical to produce 
this ronchigram. Approximately 1.5 fringes are visible. The aforementioned 
Yin-Yang pattern is visible near the center of the mirror. Figure G-3 shows the 
line spread simulation of this configuration. The wavy, horizontal shear 
suggests that some aberration will be introduced by this setup. Comparison of 
ronchigram in Figure D-l and the simulation in Figure G-3 requires some care. 
As was discussed previously in section V. B., only a very limited number of rays 
could be used to simulate the line spread function. This limitation precludes 
using spot scatter diagrams with high density, which would permit reliable line 
spread statistics to be determined. Consequently, it is noted that almost all 
calculated line spread widths are intrinsically wider than their corresponding 
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observed line spread widths. This trend is repeated for all succeeding 
ronchigrams. 

b . Figure D-2 

The entrance slit and the reticle were both rotated 30 ° from 
vertical to produce this ronchigram. The picture was overexposed, but 
approximately 1.5 fringes are visible. This indicates that the 30° rotation did 
not induce any appreciable increase in the overall system aberration. 

c. Figure D-3 

The entrance slit and the reticle were both rotated 60 ° from 
vertical to produce this ronchigram. Approximately 1.5 fringes are visible here 
also. The Yin-Yang effect is slightly visible near the center of the mirror. The 
conclusion , therefore, is that the 60° rotation does not cause a significant 
deterioration in image quality either. 

d. Figure D-4 

The entrance slit and the reticle were both rotated 90° from 
vertical (horizontal) to produce this ronchigram. Approximately one fringe is 
visible. The Yin-Yang effect is visible near the center of the mirror. Figure G-4 
shows the line spread simulation of this configuration. The wavy, vertical shear 
suggests that some additional aberration will be introduced by this 
configuration. However, note from the scale in Figure G-4 that this shear is 
less than 10% of the overall line spread width. This effect should be 
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unmeasurable with the Ronchi test used here. Indeed, no obvious evidence of 
it is seen. 

e. Figure D-5 

The entrance slit and the reticle were both rotated 120 ° from 
vertical to produce this ronchigram. Approximately one fringe is visible. This 
appears to be the best orientation for a line object because it produces the 
fewest Ronchi fringes and the least amount of aberration overall, for a single 
mirror. 

f Figure D-6 

The entrance slit and the reticle were both rotated 150 ° from 
vertical to produce this ronchigram. Approximately one fringe is visible. The 
Yin-Yang effect is visible near the center of the mirror, so this configuration is 
not quite as aberration free as the previous case. 

4. Complete System Containing 2 Mirrors and 2 Masks 

Figures E-l through E-6 show how the aberrations of the optical 
system affect the light propagating through the system with the slit and reticle 
at various orientations about the principle ray. The ronchigrams were 
photographed at the position of Mask 2. 

a. Figure E-l 

The entrance slit and the reticle were both vertical to produce 
this ronchigram. Approximately 1.5 fringes are visible. Figure G-5 shows the 
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line spread simulation of this configuration. The wavy, horizontal shear 
suggests that some additional aberration will be introduced by this setup, but 
again it is too small to detect directly with this Ronchi test. There is some 
evidence also for coma, as shown by the variable width of the spot diagram. 

6. Figure E-2 

The entrance slit and the reticle were both rotated 30 ° from 
vertical to produce this ronchigram. Approximately 1.5 fringes are visible. The 
half-moon effect is visible in this photograph. As described earlier (see section 
V.E. and Reference 13), this Ronchi pattern is indicative of both astigmatism 
and coma. 

c. Figure E-3 

The entrance slit and the reticle were both rotated 60 ° from 
vertical to produce this ronchigram. Approximately 1.5 fringes are visible. The 
half- moon effect is visible in this photograph. Note that the curvature of the 
fringes is more pronounced than in Figure E-2, indicating that even more 
astigmatism and coma are present at this orientation. 

d. Figure E-4 

The entrance slit and the reticle were both rotated 90 ° from 
vertical (horizontal) to produce this ronchigram. Approximately 1.5 fringes are 
visible. The most distorted half-moon pattern is visible in this photograph, 
implying the most severe coma and astigmatism are seen at this orientation. 
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Figure G-6 shows the line spread simulation of this configuration. The wavy, 
vertical shear suggests that some aberration will be introduced by this setup. 
More obviously, the left margin of the calculated line spread function is 
substantially wider than is the right margin. This is characteristic of 
substantial coma, and is entirely consistent with the experimental observations. 
Figure G-7 shows the line spread simulation in Figure G-6 using 1100 uniformly 
distributed rays to illuminate the mirror. The other line spread functions in 
Appendix G also have a similar shape when they are illuminated in this manner. 

e. Figure E-5 

The entrance slit and the reticle were both rotated 120 ° from 
vertical to produce this ronchigram. Approximately 1.5 fringes are visible. The 
half-moon effect is not as extreme in this photograph as in the previous ones. 
Similar to the ronchigram in Figure D-5 (produced at the same orientation), this 
ronchigram shows the least aberration and the best focus of the ronchigrams. 

f. Figure E-6 

The entrance slit and the reticle were both rotated 150 0 from 
vertical to produce this ronchigram. Approximately 1.5 fringes are visible. The 
half-moon effect is visible in this photograph, at about the same level as in 
Figure E-5. 
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F. FINAL ANALYSIS 



It is therefore concluded that the system aberrations, as predicted by the 
calculations and measured by the Ronchi test, are within the acceptable 
tolerances. The basic optical design of the multiplexing imager can successfully 
incorporate simple spherical mirrors, without introducing a noticeable 
degradation in system performance. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



As stated in the introduction, the primary goal of this thesis was to 
address the following questions: (1) Can a small object be reimaged accurately 
onto two-sided encoding Walsh masks using inexpensive, concave, spherical 
mirrors? (2) If so, what is the optimum configuration to achieve this? These 
questions have been answered successfully. The prototype optical system used 
two mirrors with radii of curvature = 128.70 cm and R 2 = 152.30 cm and 
their locations specified in Table 4. When the system dimensions are scaled, the 
prototype was sightly larger than the original goal for the imaging system of one 
square meter. Once the actual disks with the transmitting-reflecting encoding 
Walsh masks, the chopper and the positioning system are developed, the spacing 
between the optical elements can be measured and used to refine the optical 
system. This may well reduce the radii of curvature needed for the mirrors 
used in reimaging, thus reducing the overall dimensions of the imaging system. 

The initial design specification to keep the RMS point spread function less 
than 50 pm was met. This kept the aberrations in the optical system to a 
minimum. In a previous thesis by LT Musselman [Ref. 3], it was shown that 
diffraction should not seriously effect the multiplexed imaging instrument at the 
middle-infrared wavelength. This thesis research complemented that work, 
showing that geometric aberrations will not preclude the development of such 
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an instrument. Although small amount of astigmatism and coma were evident 
in the ronchigrams photographed at the positions of mask 1 and mask 2 (see 
Appendices D and E), it may be possible to reduce these aberrations if 
aspherical mirrors are used. However, it does not appear that the aberrations 
will significantly degrade the instrumental performance. 

It is recommended that further thesis research be conducted to continue 
the full development of an infrared imaging system. 
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APPENDIX B 




Figure B-l Ronchigram of Mirror 1 at a Distance Equal to its Radius of 
Curvature (R^ from the Reticle and the Entrance Slit 
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Figure B-2 Ronchigram of Mirror 1 at a Distance Equal to its Radius of 
Curvature (Rj) + 5 mm from the Reticle and the Entrance Slit 
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Figure B-3 Ronchigram of Mirror 1 at a Distance Equal to its Radius of 
Curvature (R t ) + 10 mm form the Reticle and the Entrance Slit 
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Figure B-4 Ronchigram of Mirror 1 at a Distance Equal to its Radius of 
Curvature (Rj) + 15 mm from the Reticle and the Entrance Slit 
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APPENDIX C 




Figure C-l Ronchigrams of Mirror 2 at a Distance Equal to its Radius of 
Curvature (R 2 ) from the Reticle and the Entrance Slit 
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Figure C-l Ronchigram of Mirror 2 at a Distance Equal to its Radius of 
Curvature (R 2 ) + 5 mm from the Reticle and the Entrance Slit 
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Figure C-3 Ronchigram of Mirror 2 at a Distance Equal to its Radius of 
Curvature (R 2 )+ 10 mm from the Reticle and Entrance Slit 
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APPENDIX D 




Figure D-l Ronchigram at Mask 1, Vertical Reticle and Entrance Slit 





Figure D-2 Ronchigram at Mask 2, 

Reticle and Entrance Slit Rotated 30° from Vertical 
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Figure D-3 Ronchigram at Mask 1, 

Reticle and Entrance Slit Rotated 60° from Vertical 
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Figure D-3 Ronchigram at Mask 1, 

Reticle and Entrance Slit Rotated 90° from Vertical 
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Figure D-5 Ronchigram at Mask 1 
Reticle and Entrance Slit Rotated 120° from Vertical 
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Figure D-6 Ronchigram at Mask 1, 

Reticle and Entrance Slit Rotated 150° from Vertical 
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APPENDIX E 




Figure E-l Ronchigram at Mask 2, Vertical Reticle and Entrance Slit 





Figure E-2 Ronchigram at Mask 2, 

Reticle and Entrance Slit Rotated 30° from Vertical 
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Figure E-3 Ronchigram at Mask 2, 

Reticle and Entrance Slit Rotated 60° from Vertical 



71 





Figure E-4 Ronchigram at Mask 2, 

Reticle and Entrance Slit Rotated 90° from Vertical 
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Figure E-5 Ronchigram at Mask 2, 

Reticle and Entrance Slit Rotated 120° from Vertical 
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Figure E-6 Ronchigram at Mask 2, 

Reticle and Entrance Slit Rotated 150° from Vertical 
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APPENDIX F 



KEY FOR APPENDIX F 

SPOT1.XXX tables were used to produce line spread functions to 
determine the radius of curvature for Mirror 1 (R = 60”). SPOT2.XXX tables 
were used to produce line spread functions to determine the radius of curvature 
for Mirror 2 (R = 50”). The SPOTS.OPT table was used with SPOTS.RAY 
(horizontal line) and SPOTZ.RAY (vertical line) tables to produce a line spread 
function at Mask 2. SPOTS l.RAY (horizontal line) and SPOTZ1.RAY (vertical 
line) tables were used to produce line spread functions at Mask 1. Ray tables 
with an "a" (XXXXXa.RAY) were used for P45, ray tables with a "b” 
(XXXXXb.RAY) were used for P25i, and ray tables without an "a" or "b” were 
used for P25o. 
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Figure G-l Simulated Line Spread Function for P45 Rj Mirror 1 
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Figure G-2 Simulated Line Spread Function for P45 Mirror 2 



90 



0 . 2 _ 



0 . 1 . 



0 . 0 . 



- 0 . 1 - 



- 0 . 2 . 



- 0.02 - 0.01 
yf 

SP0TZ1A . RAY 



Figure G-3 Simulated Line Spread Function for P45 Mask 1 Vertical Slit 
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Figure G-4 Simulated Line Spread Function for P45 Mask 1 Horizontal Slit 
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Figure G-5 Simulated Line Spread Function for P45 Mask 2 Vertical Slit 
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Figure G-6 Simulated Line Spread Function for P45 Mask 2 Horizontal Slit 
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Figure G7 Simulated Line Spread Function for P45 Mask 2 Horizontal Slit 
1100 Uniformly, Randomly Distributed Rays 
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